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The deformation of Ne isotopes in the island-of-inversion region is determined by the double- 
folding model with the Melbourne p-matrix and the density calculated by the antisymmetrized 
molecular dynamics (AMD). The double-folding model reproduces, with no adjustable parameter, 
the measured reaction cross sections for the scattering of 28_32 Ne from 12 C at 240MeV/nucleon. The 
quadrupole deformation thus determined is around 0.4 in the island-of-inversion region and 31 Ne is 
a halo nuclei with large deformation. We propose the Woods-Saxon model with a suitably chosen 
parameterization set and the deformation given by the AMD calculation as a convenient way of 
simulating the density calculated directly by the AMD. The deformed Woods-Saxon model provides 
the density with the proper asymptotic form. The pairing effect is investigated, and the importance 
of the angular momentum projection for obtaining the large deformation in the island-of-inversion 
region is pointed out. 

PACS numbers: 21.10.Gv, 21.60.Ev, 21.60.Gx, 25.60.Dz 



I. INTRODUCTION 

Exploring the so-called "Island of inversion" is one of 
the most important current subjects in nuclear physics. 
The term "Island of inversion (IOI)" was first introduced 
by Warburton [l[ to the region of unstable nuclei from 
30 Ne to 34 Mg. In the region, the low excitation energies 
and the large B(E2) values of the first excited states 
suggest strong deformations This indicates that 

the N = 20 magic number is no longer valid. This novel 
quantum property has triggered extensive experimental 
and theoretical studies on the IOI region. 

Another important progress of research on unstable 
nuclei is the discovery of the halo structure 041] ■ Very 
recently, the interaction cross section cti was measured 
by Takechi et at for the scattering of 28 ~ 32 Ne at 
240 MeV/nucleon and it was shown that <t\ is quite large 
for 31 Ne. A halo structure of 31 Ne was reported by the ex- 
periment on the one- neutron removal reaction . This 
is the heaviest halo nucleus in the present stage suggested 
experimentally. The nucleus resides in the IOI region. 
The reaction cross section ctr (or cti) and the nucleon- 
removal cross section with radioactive beams are thus 
important experimental tools for exploring unstable nu- 
clei 0413 ■ For the scattering of unstable nuclei at in- 
termediate energies, o\ agrees with ctr exactly or nearly, 
since projectile excitations to its discrete excited states 
do not exist or small even if they exist. This is discussed 
in this paper. 

A useful theoretical tool of analyzing ctr is the mi- 
croscopic optical potential constructed by the double- 
folding model (DFM) with t he q -matrix effective nucleon- 
nucleon (NN) interaction |13l - |2lj . when the projectile 
breakup is weak. For the nucleon-nucleus scattering, the 
single- folding model with the Melbourne g- matrix well re- 
produces the data on ctr and the elastic-scattering cross 



section systematically [2(|. For 31 Ne scattering from 12 C 
at 240 MeV/nucleon, the breakup cross section is about 
1% of or j22[. The DFM is hence reliable also for other 
projectiles 20_32 Nc, since 31 Ne is a most-weakly bound 
system among them. 

In the DFM, the microscopic optical potential is con- 
structed by folding the g-matrix with projectile and tar- 
get densities. The density profile changes, if it is de- 
formed. The elongation makes the surface diffuseness and 
the root-mean-square (RMS) radius effectively larger and 
eventually enhances ctr. The amount of deformation is 
thus important. Nuclei in the IOI region are spherical or 
only weakly deformed in the Skyrme and / or Gogny HF 
(HFB) calculations; see, e.g., Refs. 0, [24[. It is even 
pointed out that the observed large B(E2; 2 + — > + ) val- 
ues can be understood as a large amplitude vibration 
around the spherical shape j25j . In such a situation, the 
additional correlations by the angular momentum pro- 
jection (AMP) often leads to possible deformed shapes; 
see Ref. [26| for Ne isotopes. 

Recently a systematic analysis was made by the anti- 
symmetrized molecular dynamics (AMD) with the Gogny 
D1S interaction for both even and odd N nuclei in the 
IOI region [H HI- The AMP- AMD calculations, i.e. 
the AMD calculation with the AMP performed, yields 
rather large deformations. This is consistent with the 
AMP-HFB calculations [Mill]- A consistent picture of 
even and odd isotopes has been obtained by the AMP- 
AMD approach, where the n-particle m-hole excitations 
of the Nilsson orbits play important roles to determine 
deformed configurations. Although it is difficult to distin- 
guish the dynamic shape-fluctuation and static deforma- 
tion in these light mass nuclei, one may use the deformed 
shape suggested by the AMD calculation to see its effect 
on ctr. Very recently the Woods-Saxon mean-field model 
with the deformation obtained by the AMP-AMD cal- 
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culation was applied for 28_32 Ne and the DFM with the 
density of the mean-field model was successful in repro- 
ducing the data on or in virtue of the deformation [29(. 

In principle, we can calculate the double-folding poten- 
tial directly by using the nucleon density calculated with 
the AMD. The nucleon density is, however, inaccurate 
in the asymptotic region, since each nucleon is described 
by a Gaussian wave packet in the AMD. Very lately we 
proposed a way of making a tail correction to the AMD 
density [3(|. Although the calculation based on the res- 
onating group method is quite time-consuming, it was 
applied to 31 Ne [3(J, since the correction is most signifi- 
cant for 31 Ne that is a most- weakly bound system among 
2 °- 32 Ne. The tail correction to a R is about 3 % for 31 Ne. 
The DFM with the tail-corrected density reproduces the 
measured or for 31 Ne, whereas the DFM without the tail 
correction underestimates the data slightly. 

In this paper, we analyze the reaction cross section 
for the scattering of 20 ~ 32 Ne from a 12 C target at 
240 MeV/nucleon by using the DFM with the Melbourne 
^-matrix in order to determine deformations of 20_32 Ne 
systematically. Here the projectile density is constructed 
either (I) by the AMP- AMD calculation with the Gogny 
D1S interaction or (II) by the Woods-Saxon model with 
the deformation obtained by the AMP- AMD calculation. 
Model I has no adjustable parameter, but the density is 
inaccurate in the asymptotic region. Model II provides 
the nucleon density with the proper asymptotic form, 
but the model includes potential parameters. As the po- 
tential parameter set, we use the parameter set recently 
proposed by R. Wyss [3l|. This set is intended to re- 
produce the spectroscopic properties of high-spin states 
from light to heavy deformed nuclei, e.g., the quadrupole 
moments and the moments of inertia, and at the same 
time the RMS radii crucial for the present analysis. 

Models I and II yield almost the same <tr for 24_29 Ne 
that have large one-neutron separation energies. Fur- 
thermore, this agreement is seen for 31 Ne, when the tail 
correction is made in Model I. This indicates that Model 
II is a handy way of simulating results of Model I with 
the tail correction. Model II is quite practical compared 
with Model I with the tail correction that requires time- 
consuming calculations. Model I with the tail correction 
and Model II reproduce the measured or of 2 °. 28 - 32 N e . 
The deformation of 28_32 Ne is then definitely determined 
through this analysis. This analysis also yields a reason- 
able prediction for the deformation of 21_27 Ne. We also 
confirm that 31 Ne is a halo nucleus with large deforma- 
tion. Furthermore, we analyze the AMP effect and the 
pairing effect on ctr by using Model II. 

The theoretical framework is presented in Sec. |TTJ We 
explain the DFM, the AMP- AMD and the Woods-Saxon 
mean-field model. We present a handy way of making 
a center-of-mass (CM) correction to the density of the 
mean-field model. We also show in Sec. |TT] that the dy- 
namical deformation effect and the reorientation effect 
neglected in the present DFM are small. This indicates 
that <7j w ctr . Numerical results are shown in Sec. IIIII 



Section [IVI is devoted to summary. 

II. THEORETICAL FRAMEWORK 

A. Double folding model 

We start with the many-body Schrodinger equation 
with the realistic NN interaction Vij for the scattering 
of projectile (P) on target (T): 

(T R + h P + h T + ]T - = , (1) 

where E is the energy of the total system, Tr is the ki- 
netic energy of relative motion between P and T, and hp 
(/it) is the internal Hamiltonian of P (T). The multiple- 
scattering theory [H, HH for nucleon-nucleus scattering 
was extended to nucleus- nucleus scattering [34j . Accord- 
ing to the theory, Eq. (fTJ) is approximated into 

(T R + h P + h T + T v - E ) ¥+) = . ( 2 ) 

iGPJGT 

where nj is the effective NN interaction in the nuclear 
medium. The Brucckncr g matrix has often been used as 
such Tij in many applications; see for example Refs. (T3 - 
l2l| . The g matrix interaction includes the nuclear- 
medium effect, but not the effect of collective excitations 
induced by the surface vibration and the rotation of finite 
nucleus, since the interaction is evaluated in the nuclear 
matter. The effect of collective excitations is small as 
shown in Sec. Ill El 

In the scattering analyzed here, furthermore, the 
projectile-breakup effect is small, since the target is light 
and E is large. This is explicitly shown in Sec. IIII Bl In 
this situation the DFM becomes reliable. In the model, 
the potential U between P and T consists of the direct 
and exchange parts, U u and £/ EX , defined by [35ll3(| 

U DR (R)=Y f p^r P )p v T (r T )^(s; Pflv )dr P dr T , (3) 

U EX (R)=Y / f Pp(rp,rp + S K(r T ,r T - S ) 

x ffifa /V) cx P [iK{R) ■ s/M}drpdr T , (4) 

where s = rp — j"t + R for the coordinate R of P from 
T. Indices \x and v stand for the ^-component of isospin 
; fx = 1/2(1/ = 1/2) means neutron and fj, = —l/2(y = 
— 1/2) does proton. The original form of C/ EX is a non- 
local function of R, but it has been localized in Eq. Q 
with the local semi-classical approximation [l5j in which 
P is assumed to propagate as a plane wave with the local 
momentum hK(R) within a short range of the NN inter- 
action, where M = AAt/(A + At) for the mass number 
A (At) of P (T). The validity of this localization is shown 
in Ref. 
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The direct and exchange parts, g^ and g^/, of the ef- 
fective NN (g-matrix) interaction are assumed to depend 
on the local density 



/V = p£{r P + a/2) + p» T {r T - a/2) 



(5) 



at the midpoint of the interacting nucleon pair. The di- 
rect and exchange parts, g^ and g®^~, are described by 

9m=±x(s; /v) = \ E S 2 g s ,l(s; P^)6^, (6) 
s 

C,t,=o(s; /v) = \Y. $ X?( s > /v)C> (7) 
-JZ^) s+1 s 2 9 s ^p^5^A^) 

s 

lj2(-l) S+T S 2 g^(s; Pflv )S^,(9) 



9^/,T z =o( s i P^v) 



S,T 



with S 1 = + 1, in terms of the spin-isospin com- 
ponents <?^J of the g matrix interaction. As for the g 
matrix interaction, we t ake a sophisticated version of the 
Melbourne interaction 
Bonn-B NN potential 



20f that is constructed from the 
38j | . In actual calculations, the 
relativistic kinematics is taken for Tr and E. 



B. AMD framework and inputs for the reaction 
calculations 

The framework and calculational procedure of the 
AMD in this study are common to those of Ref. [28| . 
and the reader is directed to it for more detail. The 
Hamiltonian of the AMD is given as 



H = T u 



(10) 



The Gogny D1S [39J is used as an effective nucleon- 
nucleon interaction Vij ; here note that the Coulomb part 
of Vij is approximated by a sum of twelve Gaussians. 
T to t and T cm represent the kinetic energies of nucleons 
and center-of-mass motion, respectively. 

The variational wave function is a parity-projected 
wave function and the intrinsic wave function is a Slater 
determinant of nucleon wave packets, 



r=PMK V2 ,..,^}: 



(11) 



where P 7 ' is the parity projector. The nucleon wave 
packet tfi is a direct product of spatial spin \i an d 
isospin £i parts, 



tfi = (j)i(r)xS 



ur)= n 



2vo 
n 



1/4 



cxp 




where the centroids of Gaussian wave packets Zi, the 



idth of Gaussian v a and the spin direction + and a 



are the parameters determined variationally as explained 
below. Note that the center-of-mass wave function is 
analytically separable from the variational wave function 
Eq. (HU: 



■i- r 



n 



sfAi 



1 A 



(15) 



(16) 
(17) 



where X represents the center-of-mass coordinate. Usu- 
ally, Z cm defined by Eq. (fT7|) can be set to zero without 
loss of generality. This is common to the angular momen- 
tum projection and GCM calculation, and all quantities 
used as inputs of reaction calculation are free from the 
spurious center-of-mass motion. 

The parameters in Eq. (jlip are determined using the 
frictional cooling method to minimize the total energy 
under the constraint on the matter quadrupole deforma- 
tion parameter /3. Here the quadrupole deformation pa- 
rameters arc defined as, 



<*») 



1/2 



K*W) <^ 2 )] 1/6 

[(a*)(y*) (z 2 )] 1/e 

vV2 



\{x 2 ) (y 2 ) (z 2 )] 1/e 



cxp 



cxp 



cxp 




5 - / 2tt 
-0COB 7+y 



5 - / 2tt 

4^ C0S (j~T 



/9cos7 



(18) 
(19) 
(20) 



Here, (x 2 ), (y 2 ) and (z 2 ) are calculated from <J>i nt in the 
intrinsic frame that is so chosen to satisfy the relation 
(a; 2 ) < (y 2 ) < (z 2 ). The constraint is imposed on the 
value of f3 from to 1.0 with the interval of 0.025. Since 
we do not make any assumption on the spatial symmetry 
of the wave function and do not impose any constraint 
on 7, it has an optimal value for each given value of /3. 

After the variation, we perform the angular momentum 
projection for each value of /3, 



r MK 



2J+ 1 
8tt 2 



dnD{; K (n)R{n), 



(21) 
(22) 



where D'{ [K (Vt) and i?(f2) are Wigner's D function and 
rotation operator, respectively. The integrals over three 
Euler angles f2 in Eq. (|2"2"|) are performed numerically. 

The AMD calculation is completed by performing the 
GCM. The wave functions that have the same parity and 
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angular momentum (J, M) arc superposed as 

<*= E E c «^) $ m^)- ( 23 ) 

In other words, K and /3 are the generator coordinates in 
this calculation. The coefficients c n K(P) are determined 
by solving the Hill- Wheeler equation. 

E H KK'0J')CnK'0') = E n J2 N K K>WJ')CnK'(P'), 
P'K' fj'K' 

(24) 
(25) 

The ground state wave function $^ ^ obtained by this 
procedure is used in the discussion of Sec. IIIII 

For the reaction calculation, two types of applications 
are performed. One is to use the deformation parameters 
P and 7 as inputs of deformed Woods-Saxon potential. 
We assign the deformation of the AMD wave function 
by picking up a GCM basis wave function that has the 
maximum overlap with the ground state wave function, 
\{$™ w \$mk(P))\ 2 - Alld then we define the deformation 
of AMD wave function as equal to P and corresponding 
7- 

The other is the direct use of the nucleon density cal- 
culated from the ground state wave function as an input 
of the double-folding potential, 



Pjmjm 



(r) = (ClE^^ X r)|4> g 7 i /V ), (26) 



2 J 



Y J P { j ) j{r){JM'X^\JM)Y^(r), (27) 



A=0 



where the summation of A in Eq. (|27[) runs for even num- 
bers. When only is taken in the double-folding po- 
tential, the resultant folding potential becomes spherical. 
This approximation is often used as a standard manner 
in the DFM. The validity of this approximation is shown 
in Sec. UTEl 



C. Woods-Saxon mean-field model 

We also perform DFM calculations with resultant den- 
sity calculated by the simple mean-field model based on 
the Woods-Saxon potential. Sophisticated AMD calcula- 
tion is a powerful tool but it costs much time to obtain re- 
liable information, so, we take a deformed Woods-Saxon 
model for alternative approach and for further investiga- 
tion. As shown later, in many cases of Ne isotope induced 
reaction cross section, the deformed Woods-Saxon model 
gives results consistent with the AMD calculation if one 
uses the consistent deformation. 



The deformed Woods-Saxon potential in Schrbdinger 
equation is composed of the central and spin-orbit part, 
which have the following forms: 



V 



V so (r) = A so 



2m rod c 



1 + cxp [distx; (r)/a] ' 



W c (r) • <r x -V 



(28) 



(29) 



where m ro d = m(A — I)/ A and dists(r) represents a dis- 
tance between a given point r and the deformed surface 
£ specified by the radius, 



R(8, 4>; a) = Roc v (ct) 



1 



A// 



0) , (30) 



with the deformation parameters a = {a\^}. The con- 
stant c v (a) is introduced to guarantee the volume conser- 
vation of nucleus. A set of deformation parameters used 
in the present work is (/?2,7, Pa) [40| . which are related 
to (a 2fl , a ifJl ) by 



£220 
«22 



Pi cos 7, 
a 2 -2 = 



75 P2 sin 7, 



V2 



"40 = |/?4(5cos 2 7_+ 1), 

"42 = «4-2 _ 
Q?44 = Ct4— 4 



6 Pi cos 7 sin 7, 
72 Pi sin 2 7, 



(31) 



where note that the other are zero. 

As for the parameter set of the Woods-Saxon potential, 
i.e., the potential depth Vo, the nuclear radius R , and 
the diffuseness parameter a of the central potential, as 
well as those for the spin-orbit potential, we take the one 
provided recently by R. Wyss [3l| (see Ref. [4l[ for the 
actual values of parameters). For proton, the Coulomb 
potential created by charge (Z — l)e that has a uniform 
distribution inside the surface £ is added to Eq. (1281) ; 
more detailed description are explained in e.g. Refs. |42j. 

The deformation parameters in the Woods-Saxon 
potential are determined by the standard Strutinsky 
(microscopic- macroscopic) method [43|, |44|, where the 
pairing correlation is included within the BCS approxi- 
mation. The monopole pairing interaction is used and its 
strength is determined according to the smoothed pair- 
ing gap method. As for the macroscopic part energy, the 
liquid-drop model of Ref. (45| is employed. The Ne iso- 
topes in the IOI region studied in the present work are 
sitting near the drip line. In such a case the standard 
Strutinsky method has problems related to the contin- 
uum single-particle states. Recently the problems are 
solved by using the so-called Kruppa prescription (46| . 
We have employed this improved method for the calcu- 
lation of both the shell correction and the pairing corre- 
lation (Kruppa-BCS method). 

As it is discussed in the following sections, we utilize 
different models for the analysis of reaction cross sections. 
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In order to compare the deformation in different models, 
it is necessary to transform the deformation parameters 
defined within each model, e.g., between (,8,7) in the 
AMD model and (,52,7, AO m the Woods-Saxon model. 
This has been done in the following way. The deformed 
surface in Eq. (f30|) defines the deformation parameter 
in the Woods-Saxon potential. We define the uniform 
density with the sharp cut surface E, 



AmiO) = p 9(R(6, </>; a) - r), 



(32) 



where po is the average density and 9{x) is a step func- 
tion, and calculate the expectation value with it, e.g., 



)uni 



•E Pnni 



(r)dr. 



(33) 



For a given AMD deformation parameters (/3, 7), the ra- 
tio of the AMD expectation values (x 2 ) : (y 2 ) : (z 2 ) is 
fixed according to Eqs. (| — (f2"U)) . Then the correspond- 
ing Woods-Saxon parameters (/?2,7, ^4) are defined to 
give the same shape, i.e, by the condition (x 2 ) : (y 2 ) : 



)uni 



{z )uni- This condition gives 



only two independent equations so that the (,#2,7) is de- 
termined under some fixed value of ,$4. We set /?4 = 
for simplicity in order to define the (/?2 , 7) values corre- 
sponding to the AMD calculation. It is then found that 
the pairs and (/?2,7) take similar values; see Ta- 

ble I VII in Sec. IIIII We have checked that the results of 
the final reaction cross sections change very little (order 
of few mb), if we use the non-zero /?4 values within the 
range -0.1 < /? 4 < 0.1. 

In the actual calculation, the Woods-Saxon potential 
is diagonalizcd with the anisotropic harmonic oscillator 
basis, where the three frequencies, uji (i = x,y,z), are 
taken to be proportional to 1 / y/ (xf) un i, which is close to 
the optimal choice. As for the basis size, we have used 
the oscillator shells A osc = n x + n y + n z = 18 in most 
cases. However, the density distribution of the nucleus 
near the drip line extends considerably, and then we have 
checked the convergence of the results carefully by taking 
larger number of shells A osc > 20 in such a case. 

As in the Hartree-Fock (HF) or Hartree-Fock- 
Bogoliubov (HFB) approach, the occurrence of defor- 
mation in the Woods-Saxon model is a symmetry- 
breaking phenomenon. The many-body wave function $ 
is then considered to be that in the intrinsic (body-fixed) 
frame [47j . and so is the nucleon density calculated with 



P 



>(r) = (*l£*(r<-r)|*) 



£ 



(34) 



as the DFM that is done in the laboratory (space-fixed) 
frame. 

One way to recover the spherical symmetry and trans- 
form the density in the intrinsic frame to that in the 
laboratory frame is to perform the angular momentum 
projection, as already explained in Eqs. (|26|) — (|27)l in the 
AMD framework. We have performed the projection cal- 
culation (without the CM correction, which is discussed 
in the next subsection) by using the method of Ref. [49| 
for the Woods-Saxon model with the BCS pairing correla- 
tion. It is found that the projected density p j°j (r) /\f^K, 

where pjj~°\r) is defined in the same way as in Eq. (|27[) 
and used in the DFM, is very similar to the following 
angle-averaged intrinsic density, 



1 

4tt 



p( in) (r, <?,</>) sin 8d8d(f>. (35) 



In Fig. [I] they are compared for 30 Ne. As it is clear, 
for both neutron and proton, the projected and angle- 
average densities are almost identical particularly in the 
tail region, while slight differences are observed in the in- 
ner region for large J values. This is rather general trend, 
and therefore, we use pSv (r) in place of p^j] (r) / \Z4tt in 
the Woods-Saxon model for the ground state of Ne iso- 
topes. 




FIG. 1: (Color online) Comparison of the projected den- 



sity, pj°j(r)/\/47r with J = 0, 2, 10, and the angle-averaged 
intrinsic density, p^v (r) in Eq. ||3SJ, for the Woods-Saxon 
model in 30 Ne. The deformation parameters are $2 = 0.4 
and 7 = /?4 = 0, and the pairing gaps A n = A p = 1 MeV. 
The results are very similar also for the case with no pairing 
correlation. 



where v ? q( t *) is the Woods-Saxon single-particle wave 
function and v a is the BCS occupation probability (the 
free contributions should be subtracted when the Kruppa 
prescription is employed, see Refs. [46|, |48[ for details). 
Therefore the deformed density p^ n '(r) = // m )(r, 9, </>) 
cannot be directly used in the reaction calculation such 



In the view of static behaviour of deformed nuclei, 
using spherical part of deformed Woods-Saxon density, 

P < j'j( r )/^/~^ ~ Pav(f)> i s we ^ justified, but, this proce- 
dure doesn't justify that dynamical aspect of deforma- 
tion is also negligible. As mentioned above, this point is 
argued in Sec. Ill El 
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D. CM correction to nucleon density of the 
Woods-Saxon mean-field model 



The projectile density is constructed with either the 
AMD or the Woods-Saxon model for Ne isotopes. In con- 
trast to the AMD calculation, the CM motion is not ex- 
cluded in the many-body wave function $ in the Woods- 
Saxon model. We then extract the CM motion from $ 
in the standard manner [50l . l5lj and propose a simple ex- 
traction prescription consistent with the standard man- 
ner. 

The wave function $ is approximated by a product of 
the CM part <E> cm and the intrinsic part $i n t: 



cm^int 



with 



A \3/4 
^b 2 ) 



A 
21P 



X' 



(36) 



(37) 



for the CM coordinate X and the size parameter b. The 
mean squared radii of $ and <j>j nt are obtained by 



<*l£ r ?l*> 



(38) 



(r 2 ) int = ($ int |^(r,-X) 2 |$ int ) (39) 



for a single-particle coordinate r^. and hence these are 
related to b as 



(r 2 ) = (r 2 ) in t 



2 A' 



(40) 



The CM correction to (r 2 ) is small (order 1/A), so it can 
be estimated with $ instead of 4>; nt : 



(^(r.-X) 2 !*) 



(41) 



The correction is a combination of the one-body and two- 
body corrections. Inserting Eq. (|4Tj) into Eq. (|40|) . we can 
determine the size parameter b and then $; n t from $. 

The proton and neutron density without and with the 
CM correction are obtained by 

p(r) = ($|£j(^-r)Pi|$), (42) 

i 

Pint(T-) = ($ int |£^(r 4 -X-r)P ? ;|$ int ), (43) 

i 

where Pi is a projector for proton or neutron. These 
densities satisfy 



p(r) = J dr'\<S> cm (r-r')\ 2 p int (r'). 



(44) 



Thus, the density pi nt with the CM correction is obtained 
by unfolding p with |$ cm | . 



Instead of the complicated unfolding procedure [501 ] . 
we can take the following simple prescription. As shown 
in Eq. (j4T)|). the difference between (r 2 ) and (T* 2 )int i s 
small, because it is of order 1/A. This indicates that r 
dependence of Pint(^) is similar to that of p{r). We can 
then approximate pi n t(r) by 



pintO) = T*p(r/a 



with a scaling factor 



= 4 1- 



3 b 2 



(45) 



(46) 



where a has been determined to reproduce (r 2 )i n t of Eq. 
(|4T|) . The error of this simple prescription to the un- 
folding procedure is only 0.1% in <tr for Ne isotopes, so 
we use the simple prescription for the density calculated 
with the mean-field model. The RMS radii, \J (r 2 ) and 
\J ( r2 )int, without and with the CM correction are esti- 
mated with the spherical HF model, and the parameter 
b is evaluated with Eq. (|40f from the RMS radii. For 
each of Ne isotopes, we use a common b among the HF 
calculation and the spherical a nd de formed WS calcula- 
tions, since the difference of \J (r 2 ) among these mean- 
field models are at most 6% and the 6% error to the 1.5% 
CM correction is negligible. 



E. Dynamical deformation and reorientation effects 

When the projectile is deformed in the intrinsic frame, 
the deformation enlarges the radius of the projectile den- 
sity in the space-fixed frame and eventually enhances the 
reaction cross section. This static deformation effect has 
already been included in the DFM by making the AMP. 
Another effect to be considered is the dynamical defor- 
mation effect that is an effect of the rotational motion of 
the deformed projectile during the scattering. This ef- 
fect on the reaction cross section is found to be small for 
intermediate-energy nucleus-nucleus scattering [29| . This 
was confirmed with the adiabatic approximation to the 
rotational motion of projectile and the eikonal approxi- 
mation to the relative motion between the projectile and 
the target. In this subsection, the effect is investigated 
with no approximation. 

In order to test the dynamical deformation ef- 
fect, we consider the scattering of 30 Ne from 12 C at 
240 McV/nucleon and do a coupled-channel calculation 
between the + ground state and the first 2 + state of 
30 Ne. The projectile density is calculated by the DWS 
model with the deformation evaluated by the AMD. The 
coupling potentials in the coupled-channel calculation are 
obtained by the so-called single-folding model; namely, 
the nucleon- 12 C potential is first evaluated by folding 
the Melbourne-^ matrix interaction with the target den- 
sity of 12 C and the coupling potentials are then obtained 
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by folding the nucleon- 12 C potential with the projectile 
transition densities. 

In the single-channel calculation with no dynamical de- 
formation effect, the resultant reaction cross section is 
1469 mb. This result overestimates the corresponding re- 
sult of the DFM by about 10 %, which is enough to accu- 
racy of the present test. In the coupled-channel calcula- 
tion with the dynamical deformation effect from the first 
2+ state, the resulting reaction cross section is 1468 mb. 
Thus the dynamical rotation effect on the reaction cross 
section is estimated as less than 0.1 %. The reason why 
the effect is small for intermediate-energy nucleus-nucleus 
scattering is shown in Ref. [29|. The integrated inelastic 
cross section to the first 2 + state is 2.9 mb. This is 0.2 % 
of ctr, indicating that <j\ ~ or. 

The folding potential U is not spherical in general, 
when the spin of the projectile is not zero. This ef- 
fect called the reorientation effect is also tested by 
the coupled-channel calculation for the scattering of 
31 Ne(3/2-) from 12 C at 240 MeV/nucleon, where the 
single-folding model is used. The resultant reaction cross 
section is 1512mb, whereas the corresponding cross sec- 
tion is 1515mb when the non-spherical part of U is 
switched off. The reorientation effect is 0.2 % and hence 
negligible for intermediate-energy nucleus-nucleus scat- 
tering. 



III. RESULTS 

A. Reaction cross sections for stable nuclei 

We first test the accuracy of the DFM with the 
Melbourne-^ matrix NN interaction for 12 C scattering at 
incident energies (Ei n ) around 240A MeV from stable tar- 
gets, 12 C, 20 Ne, 23 Na and 27 A1. Experimental data on er R 
are available for a 12 C target at E ln — 250. 8 A MeV and a 
27 A1 target at E in = 250. 7 A MeV [52|. For 20 Nc and 23 Na 
targets, a\ at E- m = 240A MeV were recently deduced 
from measured <j\ at around 1 GeV/nucleon [53|,[54| with 
the Glauber model (l0| . 

For these stable nuclei, we take the phenomenological 
proton-density [55| deduced from the electron scattering 
by unfolding the finite-size effect of the proton charge 
in the standard manner [56|, and the neutron density is 
assumed to have the same geometry as the corresponding 
proton one, since the proton RMS radius deviates from 
the neutron one only by less than 1% in the Hartree-Fock 
(HF) calculation. 

DFM calculations are done with three types of effective 
NN interactions: the Love-Franey i-matrix interaction 
(£lf) [13 1 the Melbourne- g matrix interaction (<7mp) (20| 
evaluated from the Paris realistic NN interaction [58| and 
the Melbourne interaction (<?mb) |20j constructed from 
the Bonn-B realistic NN interaction [381 ] - 

Table U shows experimental and theoretical reaction 
cross sections for a 12 C target at E ln = 250. 8 A MeV and 



TABLE I: Reaction cross sections for 12 C+ 12 C scattering 
at E in = 250.8^1 MeV and 12 C+ 27 A1 scattering at E- in = 
250. 7A MeV. Results of three types of effective nucleon- 
nucleon interactions are compared with the corresponding 
data [52|]. The cross sections are presented in units of mb. 

Target Exp. [52] t LF [57] g M p [20J g MB [20] 



2 C 782.0 ± 10.0 917.7 793.1 



795.9 



? A1 1159.0 ± 14.0 1337.5 1164.9 1185.2 



action £lf nas 110 nuclear medium effect. In this case, 
the theoretical reaction cross sections overestimate the 
mean values of data [52[ by 17% for 12 C and by 15% for 
27 A1. In the cases of <?mp and <7mb with the medium ef- 
fect, meanwhile, the overcstimation is only a few percent 
for both 12 C and 27 A1. The medium effect is thus signif- 
icant, and an amount of the effect is almost independent 
of the bare NN interaction taken. As for 27 A1, the re- 
action cross section calculated with (/mb agrees with the 
mean value of data, when the theoretical ctr is multiplied 
by the factor F = 0.978. 




700 ' — L 



27 



Al target at E in = 250. 7 A MeV. The effective inter- 



12 14 16 18 20 22 24 26 28 
A ( mass number) 



FIG. 2: (Color online) Reaction cross sections for scattering 
of 12 C on stable nuclei from A = 12 to 27. The data at 
250.8,4 MeV for 12 C and 27 A1 are taken from Ref. [H]. The 
data at 240A MeV for 20 Ne and 23 Na are deduced from mea- 
sured cti at around 1 GeV/nucleon [53l. I54T] with the Glauber 
model The solid (dotted) line stands for results of the 
DFM with gMB after (before) the normalization of F = 0.978, 
whereas the dashed line corresponds to results of £lf- 



In Fig. [2 cr R or a Y is plotted for 12 C, 20 Ne, 23 Na and 
27 A1 targets. The dotted and solid lines represent results 
of the DFM with <7mb before and after the normalization 
of F = 0.978, respectively. Before the normalization pro- 
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ccdure, the dotted line slightly overestimates the mean 
values of data for A = 20 — 27. After the normalization 
procedure, the solid line agrees with the mean values of 
data for all the targets The normalization procedure is 
thus reliable. The dashed line corresponds to results of 
£lf with no normalization. The medium effect reduces 
the theoretical reaction cross sections by about 15% for 
all the targets. 

As for the scattering of Nc isotopes on 12 C at 
240 MeV/nucleon, we perform the DFM calculation with 
<7mb and the normalization factor F. The DFM calcula- 
tion with <7mb has harmless numerical ambiguity due to 
the parameterization of <?mb; the imaginary part of the 
folding potential has a small positive value in the tail re- 
gion. If the positive imaginary part is cut, it increases the 
reaction cross section by 2 % for a 12 C projectile and by 
1 % for Ne isotopes. This cut is used in this paper. If the 
cut is not taken, F becomes 1.0 and hence the resultant 
reaction cross sections for Ne isotopes are increased by 
1 % from the present results. This numerical ambiguity 
does not change the conclusion of this paper, since the 
ambiguity is tiny. 



TABLE II: Ground-state properties of Ne isotopes predicted 
by the AMD. For 28 Ne, the oblate state with ft = 0.28 is the 
main component of the ground state, but it is strongly mixed 
by the prolate state with p — 0.5. 



nuclide J n (exp) J n (AMD) S-m [MeV] g 7 



20 Ne 


0+ 


0+ 




0.46 







21 Ne 


3/2+ 


3/2+ 


7.111 


0.44 







22 Ne 


0+ 


0+ 


9.779 


0.39 







23 Ne 


5/2+ 


5/2+ 


6.021 


0.32 







24 Ne 


0+ 


0+ 


8.231 


0.25 




60 


25 Ne 


1/2+ 


1/2+ 


4.339 


0.20 




31 


26 Ne 


0+ 


0+ 


5.153 


0.22 




0.1 


27 Ne 


(3/2+) 


3/2+ 


1.767 


0.27 


13.6° 


28 Ne 


0+ 


0+ 


3.123 


0.28(0.50) 







29 Ne 


(3/2+) 


1/2+ 


1.321 


0.43 







30 Ne 


0+ 


0+ 


2.025 


0.39 







31 Ne 




3/2~ 


0.248 


0.41 







32 Ne 


0+ 


0+ 


1.012 


0.33 








B. AMD analysis for Ne isotopes 

Table HTl shows AMD results for the ground-state prop- 
erties of Ne isotopes, i.e., the spin-parity (J ir ), the one- 
neutron separation energy S-i n and the values of defor- 
mation parameters (/?, 7). The AMD yields the same J 7r 
as the data displayed on the web site [59| , although they 
are not established experimentally for 27 > 29 > 31 Ne. Partic- 
ularly for 31 Ne in the IOI region, the ground state has 
J n = 3/2 _ and small 5_i n consistent with the data 0.290 
±1.640 [6(|. For 28 Ne corresponding to the boundary of 
the IOI region, the main component of the ground state 
is the oblate state with f3 = 0.28, but it is strongly mixed 
by the prolate state with ft = 0.5. The deformation pa- 
rameter j3 decreases as A increases from 20 to 25 and 
increases as A increases from 25 to 32. The deformation 
becomes smallest at A = 25. 

Figure [3] plots (a) the total binding energy and (b) 
S—i n as a function of A; here the data are taken from 
Refs. [H EM- In the HF and HFB calculations, the 
spherical shape is imposed with the filling approxima- 
tion, and the nucleus with A > 30 arc unbound. The 
Gogny-HF calculations (dotted lines) underestimate the 
total binding energy systematically. This situation is im- 
proved by the Gogny-HFB calculations (dashed lines). 
The Gogny-AMD calculations (solid lines) yield even bet- 
ter agreement with the data. For S-\ n , the Gogny-HF 
calculations can not reproduce the even-odd difference 
well, but this problem is improved by the Gogny-HFB 
calculations. Thus the pairing effect is important for 
S- in- The Gogny-AMD calculations almost reproduce 
the even-odd difference for all Ne isotopes from A=21 
to 32, but slightly underestimate the experimental even- 
odd difference. This may indicate that the pairing effect 



is partly included in the Gogny-AMD calculations. The 
deformation parameter /3 is 0.33 for 32 Ne and 0.41 for 
31 Ne. The reduction of /3 from 0.41 to 0.33 may come 
partly from the pairing effect. 

Figure |4] represents ctr for scattering of Ne isotopes 
on 12 C at 240 MeV/nucleon. The AMD calculations 
(solid line) succeed in reproducing the data [Io|], while 
the spherical Gogny-HF calculation (dotted line) under- 
shoots the data; note that the spherical Gogny-HFB cal- 
culation yields the same result as the spherical Gogny-HF 
calculation within the thickness of line. The nuclei with 
A > 30 are unbound in these spherical calculations. The 
enhancement from the dotted line to the solid line comes 
from the deformation of the ground state, since the defor- 
mation is a main difference between the two calculations. 
The AMD results are consistent with all the data except 
31 Ne. The underestimation of the AMD result for 31 Nc 
comes from the inaccuracy of the AMD density in its tail 
region. 

The tail problem is solved by the following resonat- 
ing group method [30]. In principle the ground state 
<&( 31 Ne; 3/2]") of 31 Ne can be expanded in terms of the 
ground and excited states $( 30 Nc; J£) of 30 Nc. This 
means that the ground state of 31 Ne is described by the 
30 Ne+n cluster model with the core ( 30 Ne) excitations. 
The cluster-model calculation can be done with the res- 
onating group method (RGM) in which the ground and 
excited states of 30 Ne are constructed by the AMD: 

$( 31 Ne;3/2]-) = £ A { Xnlj (r)Y lm (rM 30 Ne; J^ n } , 

nj-rr 

(47) 

where <f> n is the intrinsic wave function of last neutron 
and Xnii ^ s the relative wave function between the last 
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FIG. 3: (Color online) Results of the AMD, the spherical 
Gogny-HF and the spherical Gogny-HFB calculation for (a) 
the total binding energy and (b) the one-neutron separation 
energy of Ne isotopes. The dotted, dashed and solid lines 
represent results of Gogny-HF, Gogny-HFB and AMD cal- 
culations. In the spherical HF calculations, the nuclei with 
A > 30 are unbound. The experimental data are taken from 
Refs. 



FIG. 4: (Color online) Reaction cross sections for scattering 
of Ne isotopes on 12 C at 240 MeV/nucleon. The solid (dotted) 
line represents results of the AMD (spherical Gogny-HF) cal- 
culations. The dashed line with a closed square is the AMD 
calculation with the tail and breakup corrections. The ex- 
perimental data for A — 28 — 32 are taken from Ref. [To^ j - 
The data for 20 Ne is deduced from measured aj at around 
1 GeV/nucleon [k| with the Glauber model flO^ . 
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neutron and the core ( 30 Ne). Here the wave function 
of 30 Ne includes many excited states with positive- and 
negative-parity below 10 MeV in excitation energy. This 
AMD+RGM calculation is quite time consuming, but it 
was done for 31 Ne. The tail correction to <tr is 35 mb. 

For a weakly bound system such as 31 Ne, furthermore, 
the projectile breakup effect is not perfectly negligible. 
This effect is simply estimated by assuming the potential 
model for the 30 Ne+n system and solving the three-body 
dynamics of the 30 Ne+n-|- 12 C system with the method 
of continuum discretized coupled channels (CDCC) [gl, 
HH- CDCC is an accurate method for treating exclusive 
reactions such as clastic scattering and elastic-breakup 
reactions. The theoretical foundation of CDCC is shown 
in Refs. jgjjfi^. CDCC succeeded in reproducing data on 



FIG. 5: (Color online) A dependence of the coefficient C. 
Open triangles show results of the AMD calculations, whereas 
an open square corresponds to a result of the AMD+RGM 
calculation for 31 Ne. There is a non-negligible difference be- 
tween the AMD and the AMD+RGM result for 31 Ne. Thus, 
C is slightly reduced, when the density has the halo structure. 



the scattering of stable and unstable projectiles [Hll, HH, 
l67rj80T |. Here the interactions between 30 Ne and 12 C and 
between n and 12 C are constructed with the DFM with 
the Melbourne g-matrix, and the potential between 30 Ne 
and n is made with the well-depth method; see Refs. [H, 
l8l| for the potential parameters. The correction is 10 mb 
corresponding to 0.7% of ctr. In Fig. [H the dashed line 
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FIG. 6: (Color online) Matter RMS radii of Ne isotopes calcu- 
lated with the AMD, the AMD+RGM, and the spherical HF 
model. The closed circle represents results of the AMD, and 
the closed square denotes a result of the AMD+RGM model 
for 31 Ne. The opened circles are results of the spherical HF 
calculation. 



stands for the AMD result with the tail and breakup 
corrections for 31 Ne. The result is consistent with the 
data for 28 " 32 Ne. 

The r eactio n cross section or is sensitive to the RMS 
radii, yj (r 2 )p and yj (r 2 )x, of projectile and target. Ac- 
tually, the DFM calculation for 20_32 Ne projectiles shows 
that 



ctr 



Ctt 



(48) 



where C is a slowly varying function of A around C = 
12.4 ~ 13.5 [mb/fm 2 ], as shown in Fig. El Now it is 
assumed that the projectile density has a deformed well 
shape. If the volume conservation is imposed with the 
general shape in Eq. (|30[) . the matter squared radius 
(r 2 ) P of projectile is described by 



(r 2 ) P = (r 2 ) 



1 



5 

47T 



Xji 



(r 2 )o 



1 + 7^2 2 +/?4 2 + ---) 



(49) 
(50) 



up to the second order in the deformation parameters 
{ctA/i}, where (r 2 )n is the matter squared radius in the 
spherical limit. The triaxial parameter 7 does not ap- 
pear in Eq. (|50|) . This means that the triaxial defor- 
mation little affect the RMS radius and then ctr; it has 
been confirmed by utilizing the Woods-Saxon model in 
place of the AMD model, in which the triaxial parame- 
ter 7 is optimized and cannot be changed artificially. For 
27 Ne, we varied the 7 parameter from 0° to 60° with fix- 
ing P2 = 0.273 (/?4 = 0) corresponding to the /3 = 0.27 
predicted by the AMD model, but the resultant reaction 
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FIG. 7: (Color online) Neutron and proton RMS radii of 
Ne isotopes calculated with the AMD and the AMD+RGM 
model. In panel (a), the closed (opened) circle represents 
the proton (neutron) RMS radius calculated with the AMD 
and the closed (opened) square denotes a result of the 
AMD+RGM calculation for proton (neutron) of 31 Ne; note 
that the AMD+RGM result agrees with the AMD result for 
proton RMS radius. Panel (b) shows the difference between 
neutron and proton RMS radii. The open circle (square) 
stands for the AMD (AMD+RGM) result. 



cross section changes only by 0.2%. Figure [5] shows the 
RMS radii of the spherical-HF, AMD, AMD+RGM cal- 
culations. The differences among three calculations for 
the RMS radii arc similar to those for or. shown in Fig. 
H The difference between the AMD and AMD+RGM 
calculations for the RMS radius of 31 Ne is appreciable, 
indicating that the tail correction is significant for this 
very-weakly bound system. 

Finally we compare the neutron RMS radius yj (r 2 ) 

with the proton one ^/(r 2 ) in order to see the isovector 

components of the Ne-isotope densities. Figure [7] shows 

A-dependence of yj (r 2 ) and . / (r 2 ) for Ne isotopes. In 
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FIG. 8: (Color online) Neutron and proton density profiles for (a) 28 Ne, (b) 29 Ne, (c) 30 Ne, (d) 31 Ne and (e) 32 Ne on a linear 
scale and (f) 31 Ne on a logarithmic scale. In panels (a)-(f), the solid (dotted) line presents the neutron (proton) density profile 
calculated with the AMD method, whereas the dashed line in panel (f) is the neutron density profile calculated with the 
AMD+RGM method. 
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panel (a), the neutron and proton RMS radii increase 
with A, when A > 24. For A = 20 - 24, the proton RMS 
radii have a bump. This implies that at A = 20 — 22 
the proton-neutron correlation is strong and hence the 
alpha-clustering grows. Panel (b) shows the difference 



Tp) as a function of A. The difference also 

goes up as A increases. There is a sizable jump between 
A = 28 and 29, since the deformation /3 is around 0.25 
at A = 24 - 28 but around 0.4 at A = 29 - 32. As a 
result of this gap, the radius difference is around 0.35 fm, 
indicating that Ne isotopes are either skin or skin-like 
nuclei for A = 29, 30, 32. For 31 Ne, the radius difference 
calculated with the AMD+RGM method is about 0.5 fm 
that is significantly bigger than 0.35 fm. This implies 
that 31 Ne is a halo nucleus. These interpretations are 
more obvious through the neutron and proton density 
profiles shown in Fig. [51 Panels (a), (b), (c), (d), and 
(e) show the density profiles for 28_32 Ne, respectively. 
29 ' 30 ' 32 Ne have the neutron-skin structure. In panel (f), 
the density profiles for 31 Ne is plotted on a logarithmic 
scale. The neutron density (dashed line) calculated with 
the AMD+RGM method has a long-range tail, indicating 
that 31 Nc has the halo structure. 
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FIG. 9: (Color online) Reaction cross sections calculated 
for Ne isotopes calculated with the SWS model and the 
Gogny-HF method. The spherical shape is imposed for both 
the calculations. The dotted line represents results of the 
SWS model, while the solid line corresponds to the spherical 
Gogny-HF results. The nucleus with A > 30 are unbound. 
The experimental data are taken from Ref. [ic( |. 



C. Woods-Saxon mean-field model 

In this subsection, the results of the Woods-Saxon 
mean-field model are investigated. First of all, in order 
to see that the parameter set of the Woods-Saxon poten- 
tial is reasonable, the spherical case is studied. Figure [3] 
shows the reaction cross sections for Ne isotopes calcu- 
lated with the spherical Woods-Saxon (SWS) model (ne- 
glecting the pairing correlation) and the spherical Gogny- 
HF method. The SWS model (dotted line) well simu- 
lates results of the spherical Gogny-HF calculation (solid 
line). This means that the SWS model yields almost the 
same matter radius as the spherical Gogny-HF calcula- 
tion. The SWS model with the present parameter set 
proposed by R. Wyss [3l| is thus a handy way of simu- 
lating the spherical Gogny-HF calculation. 

Figure [TU] shows the neutron Nilsson diagram for 30 Ne 
calculated with the DWS model. It is emphasized that 
the relatively large shell gap with N — 20 is observed at 
the spherical shape = 0). Comparing this figure with 
the Nilsson diagram calculated with the AMD model in 
Fig. 2 of Ref. [27|, one can see that the Nilsson diagram 
of the DWS model is close to that of the AMD model. In 
both the models (AMD and DWS), the [2,0,0,1/2] and 
the [3,3,0,1/2] orbit in terms of the Nilsson asymptotic 
quantum numbers [JV,7i3,A,fi] cross each other at /3% w 
0.4, although the single-particle energy at the crossing 
point is -3 McV in the AMD model and -4 MeV in 
the DWS model. It is well-known that the occupation of 
this down-sloping orbit [3,3,0,1/2] derives the system to 
deform near N rj 20. 

Next, the results of the deformed Woods-Saxon (DWS) 



1p3/2 




FIG. 10: (Color online) The neutron Nilsson diagram for 30 Ne 
in the deformed WS model, where the other parameters are 
fixed to /?4 = 7 = 0. The solid (dashed) lines correspond to 
the positive (negative) parity orbits. The Nilsson asymptotic 
quantum numbers [N ,n3,A,Q] are attached. The number 20 
stands for a neutron magic number in the spherical limit. 



model with deformation obtained by the microscopic- 
macroscopic (Strutinsky) method are discussed. The 
pairing correlation is included within the Kruppa-BCS 
approximation, which strongly affects the resultant defor- 
mations. Two cases of the pairing strengths determined 
by the smooth pairing gap A = \2/\/~Aor 4/y/AMeV are 



13 



TABLE III: Deformation parameters (ft, ft) and the pairing 
gaps A n and A p obtained by the Strutinsky calculation with 
A = 12/VA MeV. The nucleus 31 Ne is unbound and * is 
attached. 



nuclide 


ft 


ft 


A n [MeV] A 


P [MeV] 


20 Ne 


0.000 


0.000 


2.897 


2.848 


21 Ne 


0.310 


0.048 


0.000 


2.290 


22 Ne 


0.000 


0.000 


2.801 


2.826 


23 Ne 


0.200 


0.011 


0.889 


2.596 


24 Ne 


0.000 


0.000 


2.337 


2.865 


25 Ne 


0.000 


0.000 


0.000 


2.849 


26 Ne 


0.000 


0.000 


2.157 


2.826 


27 Ne 


-0.065 


-0.002 


0.989 


2.783 


28 Ne 


0.000 


0.000 


2.345 


2.783 


29 Ne 


-0.048 


-0.001 


1.391 


2.748 


30 Ne 


0.000 


0.000 


2.286 


2.734 


31 Ne 


0.101* 


0.022* 


1.618 


2.623 


32 Ne 


0.000 


0.000 


2.416 


2.674 



considered; A = 12/y/A MeV is often used as a typical 
value for medium- heavy nuclei, while A = 4/y/A MeV is 
used to simulate the weak pairing case. The smooth pair- 
ing gap A is usually supposed to correspond to the even- 
odd mass difference. However, in the light mass nuclei, 
like the Ne isotopes considered in the present work, the 
even-odd mass differences contain considerable amount 
of the shell effect of deformed mean-field, and the weaker 
pairing correlations are suggested (HJ. The small value 
A = 4/\AA MeV is chosen in accordance with it. In 
Tables IIIII and IIV1 the results of the deformation param- 
eters (,82, ft) and the neutron and proton pairing gaps 
A n and A p are summarized; all nuclei are calculated to 
be axially symmetric in their ground states. If the stan- 
dard pairing is used, all even-even isotopes turn out to 
be spherical, which contradicts the result of the AMD 
calculation. Even with weaker pairing, 25 > 26 > 28 > 30 N e are 
calculated to be spherical; this is because of the relatively 
large N = 20 shell gap (see Fig. [TO]) , which is in contrast 
to the prediction using the tensor force in Ref. (83l |. 

In order to check the consistency of the DWS model, we 
compare the obtained deformations with the systematic 
Gogny D1S HFB calculation s 18411 . the results of which 
are available on the web site [85(. Those results are also 
axially symmetric in Ne isotopes, and the (ft, ft) de- 
formation parameters are extracted in such a way that 

(r 2 Y2o)uni/(r 2 )um and ( r4 ^o) uni/ (r 2 ) uni calculated with 
the uniform density in Eq. (|32j) reproduce the correspond- 
ing values of Gogny D1S HFB calculations tabulated in 
Ref. The resultant (ft, ft) are listed in Table El It 
can be seen that (ft, ft) calculated with the weaker pair- 
ing in Table Qv] well corresponds to those of the deformed 
Gogny D1S HFB results. This result may suggest that 
the monopole pairing strength weaker than the standard 



TABLE IV: The same as TablelUbut with A = 4/v^4 MeV. 



nuclide 


ft 


ft 


An [MeV] A 


P [MeV] 


20 Ne 


0.336 


0.111 


0.000 


0.000 


21 Ne 


0.349 


0.079 


0.000 


0.000 


22 Ne 


0.362 


0.051 


0.000 


0.000 


23 Ne 


0.291 


0.052 


0.000 


0.000 


24 Ne 


0.186 


0.018 


0.000 


1.434 


25 Ne 


0.000 


0.000 


0.000 


1.657 


Of, H T 

26 Ne 


0.000 


0.000 


0.000 


1.654 


27 Ne 


0.122 


-0.001 


0.000 


1.471 


28 Ne 


0.001 


0.000 


0.966 


1.638 


29 Ne 


0.067 


-0.001 


0.000 


1.572 


30 Ne 


0.000 


0.000 


0.000 


1.611 


31 Ne 


0.291* 


0.106* 


0.000 


0.000 


32 Ne 


0.278 


0.098 


0.922 


0.000 



TABLE V: Deformation parameters ($2, ft) extracted from 
the results of deformed Gogny-DlS HFB calculations of 
Refs. [84], [1^]. The nucleus 31 Ne is unbound and no data 
are available. 



nuclide 


ft 


ft 


20 Ne 


0.325 


0.108 


21 Ne 


0.370 


0.085 


22 Ne 


0.355 


0.016 


23 Ne 


0.234 


0.011 


24 Ne 


0.179 


0.011 


25 Ne 


-0.047 


0.001 


26 Ne 


-0.002 


0.000 


27 Ne 


-0.073 


-0.005 


28 Ne 


-0.006 


0.000 


29 Ne 


-0.060 


-0.003 


30 Ne 


-0.002 


0.000 


31 Ne 






32 Ne 


0.246 


0.096 



one is more appropriate in light nuclei in accordance with 
the analysis in Ref. [Hf . 

The calculated binding energies and one-neutron sep- 
aration energies in the DWS model are compared with 
experimental data in Fig. 1111 The binding energies 
are slightly overestimated for unstable isotopes 27_30 Ne, 
which is mainly due to the liquid drop energy (45j em- 
ployed in the present work. As for the one-neutron sepa- 
ration energies, S-i n , the calculations with the standard 
pairing strength with A = 12/vC5 MeV nicely reproduce 
the experimental data, although the weaker pairing with 
A = 4/vG4 MeV may yield slightly better agreement for 
30_32 Nc. Similar agreement is seen in the AMD, but the 
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FIG. 11: (Color online) Results of the DWS model including 
the pairing correlation with A = 12/y/A or 4/y/A MeV for (a) 
the total binding energy and (b) the one-neutron separation 
energy S-i n of Ne isotopes. The solid line is a result of A = 
12/VA MeV and the dashed line of A = 4/ -/A MeV. 



predicted deformations of Ne isotopes in the two mod- 
els are very different; those in the AMD are larger par- 
ticularly in unstable isotopes 28_32 Ne. This is because 
the binding energy reflects many structural effects. The 
binding energy is thus not a good indicator of deforma- 
tion. 

Of course, the well-known indicator of nuclear defor- 
mation is the rotational spectra and the E2 transition 
probabilities. They are, however, difficult to measure in 
unstable nuclei. As already discussed in the previous 
subsection, the reaction cross section ctr can be utilized 
instead. Figure [T2l shows or for Ne isotopes calculated 
with the DWS model. The deformation parameters and 
pairing gaps are calculated either with A = 12/y/A or 
4/v / A MeV, and listed in Table |ffl] or [TV] A result of the 
DWS model with the deformations calculated with the 
deformed Gogny-DlS HFB method HUl (Table E} is 



FIG. 12: (Color online) The reaction cross sections for Ne 
isotopes calculated with the DWS model with some sets of 
deformation parameters. The solid and dashed lines repre- 
sent results of the Strutinsky method with A = 12/y/A and 
4/\/A MeV. In the dotted line, the deformation is determined 
by the deformed Gogny-DlS HFB calculation, for which the 
pairing correlation is neglected. The experimental data are 
taken from Ref. [Tol |. 



also included as a dotted line. Comparing with the defor- 
mation parameters obtained with the AMD (Table [TTJ) , 
the used value of deformation parameters are zero or 
small around 30 Ne since the neutron number N = 20 
is a spherical magic number. As a consequence of this, 
the reaction cross sections based on both the microscopic- 
macroscopic (Strutinsky) and the deformed Gogny-DlS 
HFB methods underestimate the experimental data for 
28-31 Ne. Note that the effect of finite pairing gap alone 
is small if deformation is fixed as it is shown in the next 
subsection (see Fig. 111))) . 

The reason why the AMD calculation gives large de- 
formations in the IOI region is that the optimum defor- 
mation is searched after the angular momentum projec- 
tion (note that the same Gogny PIS force is used in the 
HFB calculation of Refs. [8J, |85j). It is known that the 
potential energy surface as a function of quadrupole de- 
formation is rather shallow for nuclei in the IOI region. 
In such a case, the energy gain of the AMP at large de- 
formation can easily change the equilibrium deformation, 
see e.g. Ref. [26| for Ne isotopes. The angular momen- 
tum projection is thus important for the IOI region to 
obtain large deformations. 

D. Woods-Saxon model with AMD deformation 

In the previous subsection, it has been shown that the 
microscopic-macroscopic method with the DWS model as 
well as the deformed Gogny-DlS HFB approach (with- 
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TABLE VI: Deformation parameter /?2 and 7 deduced from 
the AMD intrinsic density. Those with higher multipoles A > 
2 are not included. The Nilsson asymptotic quantum numbers 
of last neutron are included as the last column for axially 
symmetric cases. 



nuclide 


P 


7 




7 


[N ,ri3,A,£i\ for last-n 


20 Ne 


0.46 


0° 


0.479 


0° 


[2,2,0,1/2] 


21 Ne 


0.44 







0.456 







[2,1,1,3/2] 


22 Ne 


0.39 







0.400 







[2,1,1,3/2] 


23 Ne 


0.32 







0.325 







[2,0,2,5/2] 


24 Ne 


0.25 




60 


0.258 




60 


[2,0,0,1/2] 


25 Ne 


0.20 




31 


0.202 


31.5° 




26 Ne 


0.22 




0.1 


0.221 







[2,1,1,1/2] 


27 Ne 


0.27 


13.6° 


0.273 


14.1° 




28 Ne 


0.50 







0.526 


O 




[3,3,0,1/2] 




0.28 


60 


0.291 


60° 


[2,1,1,3/2] 


29 Ne 


0.43 







0.445 







[2,0,0,1/2] 


30 Ne 


0.39 


O 




0.400 


O 




[2,0,0,1/2] 


31 Ne 


0.41 







0.422 







[3,2,1,3/2] 


32 Ne 


0.33 







0.335 







[2,0,2,3/2] 



out the AMP) do not give expected large deformations in 
the IOI region. Therefore, we employ, throughout in this 
subsection, the deformations obtained by the AMD in 
the DWS model and compare the results with the AMD 
calculations and experimental data. In most of this sub- 
section, the pairing correlation is neglected; its effect is 
discussed at the end. 
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FIG. 13: (Color online) Matter RMS radii of Ne isotopes 
calculated with the DWS, AMD and AMD+RGM models. 
The opened circles represent results of the AMD model, and 
the closed square denotes a result of the AMD+RGM model 
for 31 Ne. The opened inverted triangles are results of the 
deformed Woods-Saxon calculation. 



Table IVII lists up the deformation parameters fa and 
7 deduced from the corresponding AMD values f3 and 7 
and used in the following DWS calculations. In Fig. HH 
the matter RMS radius calculated with the DWS model 
is compared with those of the AMD and AMD+RGM cal- 
culations. The DWS model well reproduces the matter 
RMS radii of the AMD calculation for 24 ~ 29 Ne in which 
S- i n is large. For 30_32 Nc in which S-\ n is small, the 
matter RMS radii of the AMD calculation are slightly 
smaller than those of the DWS model. The deviation 
may come from the fact that the AMD density is inac- 
curate in its tail region, since the DWS model almost 
reproduces the matter RMS radius of the AMD+RGM 
calculation for 31 Ne. For 20 ~ 23 Ne, the matter RMS radii 
of the AMD calculation are larger than those of the DWS 
model. This may imply that the a clustering is well de- 
veloped in the AMD calculation and the different type of 
deformation from those included in the present Woods- 
Saxon model, e.g., the octupole deformation («3 M ), may 
be important for 20_23 Ne. 

The nucleon density distributions are plotted in Fig. 
Hilfor 24 Ne and 31 Ne. The AMD densities (dotted curves 
decrease with increasing r more rapidly than the densi- 
ties (dashed curves) of the DWS model. The deviation 
between the two densities at large r is rather small for 
24 Ne where S-\ n is large. The deviation is, however, 
enlarged for 31 Ne in which S-\ n is small. The AMD den- 
sity is thus inaccurate at large r particularly for 31 Ne. 
A tail-correction to the AMD density can be made by 
the AMD+RGM calculation. The density (solid curve) 
has actually a long-range tail and consequently becomes 
close to that of the DWS model. As an important result, 
the density of the DWS model almost agrees with that 
of the AMD+RGM calculation for 31 Ne with fa = 0.422. 
This result indicates that the DWS model with the AMD 
deformation is considered to be a handy way of making 
a tail correction to the AMD density. 

In Fig. [151 the reaction cross sections are calculated 
for Ne isotopes with the DWS, AMD and AMD+RGM 
models. The differences among the three calculations are 
similar to the corresponding differences for the matter 
RMS radii shown in Fig. [131 as expected. Note that 
the result of the DWS model nearly agrees with that 
of the AMD+RGM model. As an important result, the 
reaction cross sections calculated with the DWS model 
are consistent with the experimental data [lol ]. 

Finally, the effect of pairing correlation is investigated 
for the deformations given by the AMD model. Two 
cases of the pairing strengths given by the standard value 
A = YljyfA MeV and the weaker value A = 4/-/A MeV 
are considered as in the previous subsection, and the re- 
sult is shown in Fig. [T5] With the AMD deformations 
the even-odd effect on the separation energies S-\ n for 
29_32 Ne is too much enhanced if the standard pairing 
strength is used, while a better fitting is obtained with 
the weaker pairing strength, as is shown in Fig. I16f a) . In 
Fig. fTBT b). the reaction cross sections for Ne isotopes are 
evaluated by the DWS calculations with different pairing 
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FIG. 15: (Color online) Reaction cross sections for Ne iso- 
topes calculated with the DWS, AMD and AMD+RGM mod- 
els. The dotted line represents results of the DWS model, 
while the solid line corresponds to the AMD results. The 
closed square represents a result of the AMD+RGM calcu- 
lation without breakup contribution. The experimental data 
are taken from Ref. [1(3 ]. 



really occurs and reflects the reaction cross section. 

In the DWS model that corresponds to the AMD model 
with the tail correction, ctr, for 32 Ne is slightly larger than 
that for 31 Ne as shown in Fig.[15j but the RMS radius for 
32 Ne is smaller than that for 31 Ne as presented in Fig. [13l 
The reduction of the RMS radius comes from that of /3a- 
It is interesting to consider what causes the reduction of 
p2- This is another interesting future subject related to 
the "pairing anti-halo" effect mentioned above. 



FIG. 14: (Color online) Density distributions of (a) 31 Ne and 
(b) 24 Ne. The dotted line represents a result of the AMD 
model, whereas the dashed line corresponds to a result of the 
DWS model. The solid line is a result of the AMD+RGM 
model. 



strengths. The reaction cross sections are enhanced a bit 
by the pairing effect, but its effect is small even with the 
standard pairing strength. 

For weakly bound systems, it is speculated that the 
pairing correlation leads to an extra binding of halo orbit 
and makes the nuclear radius shrink; it is called the "pair- 
ing anti-halo" effect [HI, |87[- Our Kruppa-BCS method 
can produce the anti-halo like effect [48| , but a reduction 
due to the pairing effect is not observed in the present 
calculations. Possible reasons may be the large deforma- 
tions, which tend to prevent the anti-halo effect, and that 
the binding of the last neutron orbit is not weak enough. 
As a future work, the deformed HFB calculation is highly 
expected to answer whether the "pairing anti-halo" effect 



IV. SUMMARY 

We determined deformations of 20_32 Nc with the fully- 
microscopic AMD model that has no adjustable parame- 
ter. The quadrupolc deformation parameter determined 
is around 0.4 in the IOI region and 31 Ne is then a halo 
nuclei with large deformation. 

We have also performed the microscopic-macroscopic 
(Strutinsky) calculations with the Woods-Saxon poten- 
tial, and found that the obtained deformations are too 
small, which is consistent with the deformed HFB calcu- 
lations without the AMP in this region. 

As a reaction model, we used the DFM with the Mel- 
bourne g-matrix. The microscopic reaction model yields 
good agreement with the measured or for 2 °. 28 - 32 N e , 
if the projectile density is constructed either (I) by the 
AMP- AMD calculation with the Gogny D1S interac- 
tion or (II) by the Woods-Saxon mean-field model with 
the deformation obtained by the AMP- AMD calculation. 
Method I has no adjustable parameter, but the tail of 
the density is inaccurate. We then made a tail correction 
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to the AMD density for 31 Nc by using the AMD+RGM 
method. Method II provides the nucleon density with 
the proper asymptotic form, but a parameter set of the 
Woods-Saxon potential should be carefully chosen. The 
parameter set recently proposed by R. Wyss [3l| is shown 
to be very successful: Method II yields almost the same 
ctr as the AMD+RGM method. This means that Method 
II is a handy way of simulating results of Method I with 
the tail correction. 



The two types of DFM well reproduce the measured 
ctr for 20,28-32^ Deformations of 28_32 Ne were defi- 
nitely determined through this analysis. This analysis 
also yields a reasonable prediction for deformations of 
21 ~ 27 Ne. We also showed that the AMP is essential to 
obtain the large deformations required for reproducing 
the measured ctr for 28_32 Ne but that the effect of the 
BCS-type pairing correlation is small there. 
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FIG. 16: (Color online) The pairing effect on (a) the one- 
neutron separation energy S-in and (b) the reaction cross 
sections or for Ne isotopes. The solid line is a result of the 
DWS model without pairing correction and the dashed (dot- 
ted) lines stand for a result of the DWS model including the 
pairing correlation with A = 4/v / ~4 (12/y/A) MeV. The ex- 
perimental data are taken from Ref. [6(| [6lJ in panel (a) and 
from Ref. [l(| in panel (b). 
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